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ON THE MOTION OF HYPERION. 

By Prof. Ormond Stone, University of Virginia, Va. 

I. If we assume that the planes of the orbits of Titan and Hyperion 
coincide, the differential equations for the motion of the latter may be written 

2 dw 
dt 



d} 

dl 



= C+m'^ fsrdf, (i) 

>- Cdzv') ^ , ^ /72D / \ 

''~''['dFj +p = ^'^^; (2) 



in which r and w are the radius vector and longitude in orbit of Hyperion, 
k the constant of the Saturnian system as derived from the orbit of Titan, 
m' the mass of Titan, in'f^R and m't^S the components, in the plane of the 
orbit, of the disturbing force in the direction of and perpendicular to the radius 
vector, and C is a constant of integration. 

We shall assume that r, w, R, and 5' can be developed into series of sines 
and cosines of multiples oi d = I' — / — X, having coefficients which are periodic 
functions of -q=- aJ- — 3/', in which / and /' are the mean longitudes in orbit of 
Hyperion and Titan, and ^ is a constant to be derived from observation, tj is 
evidently the longitude of the mean conjunction point of the two satellites, and 
/' — / their mean elongation. 

If we put 

k 



and 



kVa{\—v)^C^ m'pJ(Sr\dt; 



in which a, the mean distance of Hyperion from Saturn, is a constant, and (Sr\ 
is that portion of Sr which is independent of 6. v is thus a function of ;f only. 
We may also put 

drj = (4« — 3«') dt, 



and C=^^Va{i — v„); 

whence, if a', the semimajor axis of Titan's orbit, be taken as the unit of 
distance, we may write 

V/(I - v) = ,/(i - v„) + --yjf~~,. "\ f-T fi^*-^o^1' 

\/{i + m') 4« — 3«'|/a J '" ' 

which may be solved by a simple integration. 
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We may now write 

r^-,- =k\^a{\ —i')-\ ; 7-7 U, (3) 

di ^ ' ' I -{- m' n' — n 

in which U is that portion of / SrM which remains after omitting those terms 

which are independent of d. As is evident, 

dd — {n' — n)dt. 

Substituting in (2) the value oi dwjdt derived from (3), the former equation 
becomes 

5-^^ + ^ = ./^/', (4) 

dr r' r^ a' 

where, neglecting terms in m'^, 

P^ a-R + ^4i^ii:) ^— [/, (S) 

in which 

r 

a 



If we multiply by a^^, (4) becomes 

dh , s — I + V 
in which 



as , s — I + V in {fi\ 



in'— n)l 

2. As the equations of motion of an intermediate orbit, we shall assume 

r, = a,s,, 

3 («'-«) =—3' 

i. e. a Keplerian ellipse with mean motion 3(«' — «), semimajor axis a,, and 
eccentricity e; a, is a constant, and e a function of 7] only, satisfying the condi- 
tion that the coefficient of cos 3^ in 8s is zero. 
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Neglecting powers of e greater than the fourth, we may write 

•f / = I + i ^ 

+ e(\ — -|^) COS 3^ 

— J^(i — 1^2) cos6^ (9) 
+ f ^ co.s 9^ 

— \e^ COS 12^, 
and 

>^'(' 7^) = I - 2e{\ - \^) COS 3^ 

+ f^'(i-H^)cos6^ 

— Jj^ ^ cos 9^ 

+ \'!^^<COS 12^. 

Subtracting (8) from (6), neglecting powers and products of ^i = J — .^r and 
^v =: V — ^^ and putting 

we have 

= m'Ii{Pi cos 2^ + Qi sin 2^), 



in which 8fi = /i — ^, ^v satisfies the condition that ds shall contain no term 
independent of ^, i is any positive whole number, and P^ and Qi are assumed to 
be functions of -y only. 

3- Considering, for the present, coefficients of cos od, cos ^0, cos 6d, and 
cos gd only, we may write 

/? = I + 3^ 

— 6^(i + f^) cos 3^ 

+ 1 2^^ cos 6^ 

— ^e' cos 9^. 
If we assume 

00 

ds = 2'i(flj cos id + ^i sin id), (ii) 
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and substitute for s, and ds their values given by (9) and (ji), remembering that 
rtj = o, we have 

-5-2 = — 36^8 cos 66 — 8 1 a, cos gd, 

^1 = _ 9^(1 _ 3 ^) cos id + lie" cos 66, 
^ds = 6e\ 

— [3^(f + ¥0«6 — 6^«9] cos s6 
+ [(i + 3^^)«6 — ^^^9] cos 66 

— (3^<^6 — «») cos g6, 

— 3^(1 +1^") cos 26 
+ f ^ COS 6(?. 

Substituting in (10) and separating by indeterminate coefficients, we obtain 

- 9^(1 - |^)<?/i - 3^(1 + ¥^)«6 + 6f^«9 - 3''(l + f O^^ =- '^'^a, 

i8/5/i — 3 (I — / + 1 2dfi)ag — ya, + f ^% = m'P^, 

— ya^ — 8«5 = m'P^; 

or 

^V = ( I — f ^') i'«'/'o — 6r'«6 . 

gedix = - 3^(1 + ie')dv - (i + |e^);«'/'s - 3''(' + W^;«6 + 6f'«9, 
«5 = I Mv + 6/5/i — |( I + ^2 — 1 2^//) m'P^ — m„ 

whence, solving by trial, 

m'-'8v = (i — f ^)/'o + s^-^e. 
9;«'-l^;,.,= _ 3^(1 _i^)/^^_(i - J^^)P3 + ^(i -f.^ + i2dii)P,-ie'P„ 

m'-'a, ■= - i^P„ - ieP, _ ^ (, - |^ - i2^/z)/'e + i^/',, ^'^^ 

fn'-\= lePf — ^P^. 
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For the succeeding terms, we may write 



I —I' 



(13) 



4. Considering cosine terms which involve only those multiples of d which 

are not also multiples of 3^, equation (10) may be written, with sufficient 

accuracy, 

1 ^ff + ^bs = m'SiPi cos id, 
au 

in which we may assume 

8s =■ «i cos 6 -\- a^ cos 2d + a^ cos 4^ + . . . 
and 

/? = I + 2^83 cos 3^ + 2j8e cos 60 ; 

whence, separating by indeterminate coefficients, we have 

f «1 + /?3«2 + /?3«4 + /?6«6 + /?6«7 = »i'Pl< 

/?3«1 + ■|«2 + ft«4 + /?3«5 + i^e^S = '«''^2. 

M + /?6«2 — ^«4 + M + i^e^lO = '«''^4. 

/?6«1 + i93«2 — ¥-«6 + /?3«8 + ^6«U = '^'^S. 



)56«l + M — ¥«7 + Mo + i'e^w == ^'^7 



?• 



/36«2 + M — ¥«8 + Ml + Mu = »i'Pi, 



(H) 



Puttin pi = in'Pi, and calling J the determinant whose constituents are 
the coefficients of a,, a^, etc.; J, the determinant formed by substituting /„ p^, 
etc. for the first column in A, i. e. for the coefficients of a^ ; A^ the determinant 
formed by substituting /i, A. etc. for the coefficients of «2. etc., we have 



^1 



«4 ^ -^ ; etc. 



To expand these determinants, let us consider the determinant in its most 
general form. If we put 



^yfc.. 



I 



Uji 



a» a. 



'■jit 
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in which i </ < ^ < . . . , we have 

«51 «91 



47 



*11 "21 "31 



^12 ^22 ^32 



'J^IS «23 ^J^ss 



U-ii ^90 ^o( 



(I + lA^ + 2'^(,, + SA^,, + ...), (is) 



as may readily be seen by transforming and generalizing the method for 
expanding in products of the leading constituents given in Burnside and Pan- 
ton's Theory of Equations, p. 246. 
In the expansion of d we find 

SA^ = -A^,'+ ..., 



in which 






A= l-f-f.A + A-A+--- 




-f-f + I-A + A-A + --- 




QO 

= i'.. 



r ' '1 


I I 


li-(i-ifi-(t + ^y I- 


-{i+iy'i-ii + ifj 




r 3 


' 1 3 


=*f 


1 3 .. I 




' I + («• + *) I - U + i) 


I + (2 + 4) ' i-U + i) 



= 0; 

whence, with sufficient accuracy. 



= |y7,[i-(.-i/][i-(r+i)^]. 

Expanding .^,, we obtain 



^1 =/,•!• -*•-¥• 



A /i 

A A J 
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whence 

m'-^a,= f[i+fa-^)/?3^/'. 

in which the coefficient of ^^ is equal to A^ after the omission of those terms in 
A^ which contain as factors the reciprocals of such constituents as are situated 
at the intersection of the /-column and the diagonal; in other words, since 
A,^ = o, it is equal to the sum of the terms mentioned taken with opposite signs. 
In general, when i is not a multiple of 3, we have 



9 



-^—fi + -9- f ^ I 9 1 02~|p 

- i-^l ^ g-Plg- (, _ 3)^ + 9 - (^ + 3f J ^' J ' 

_9 9 » p. . _ 9 9 



/?3/'.-3 - ;r^.---7f^^. ^A+. (16) 



9_/2-9_(^_3)2'5 -3 9_,2-g_(,^3,2 



Psfi-6 — Ta'- /■ I ^^2 A^-+6. 



in which 



A =-3^(i+f^^), 

/?6 = 6^. 

5. For the terms involving only sines of multiples of 3^, equation (10) may 
be written 

F -jfj2 + fi^^ = ^'^'iQsi sin std, 

in which we may assume 

<^-^ = ^'i ga sin Zi6, 
1 

^ = /?o + 2li^ii COS lid, 

whence, separating by indeterminate coefficients, we have 

hgi i- (A — ?i )gi -h 0?6 — A2)<?-9 + (/?9 — ;^i5),?-i2 + • ■ • = ^3. 

{?% - ^9 ).?-3 + Kg, + (/Ss - A.0.f9 + (A - ^9)^12 + ... = ^76- 
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(^6 - ^uUs + (/?3 - A5)^6 + ^^9 + (/?3 - Al)^12 + • • • = ?9. 

0^9 — fi.d^Z + (/?6 — fiis)^e + (/?3 — Al).f9 + ''iS.^-lJ + • • • = ?12, 

where, if (T ^ 9^/z, 

^3 = - 3^ - <^ + ¥ ^S 

K = — 3 + 3^ — 4<^. 
^ =-8, 

''12= — 15- 
If D be the determinant whose constituents are the coefficients, of ^3, g-^, 
etc. ; Z^j the determinant found by substituting ^3, q^, etc. for the first column in 
D, i. e. for the coefficients of ^3, etc., we have, as in the preceding section, 

g^=-^; g.= ^; etc. 
Applying formula (15), we obtain 

D = ;3U-.- f I - ^^3-^9)' _^ _ K-ll + ?m + _Jl_ 

K, /jXj /3/g /jXg /g;^j2 ^3/6/9 ^«6^9^12 

A = M^ ...fi — ^« — ^»-:^ — -^ + ?^0 , 

V ^6?3 ^9?3 '6^9 ^9^12 ^6'9?3 J 

A =/i3M.-.ri-w. 

D —llal r I _ ^3' _ /^«^3 I 2/?3V3l 
U^ ''3''6y9''l2 • • • I YT 1 ^ Tl ' 

A2 = W9?i2^i5 • • • [ " ~ n" J ' 



whence, substituting for .^3, ^^3, etc. within the parentheses their values in terms 
of a and e, 

A = K\--- [(i — -A^) ?3 — eq^ + |?9l, 

A =^ V,2 . . . [3^^3 — (3^ + <^) ^'sj. 

A = ''6''l2''l5 . . . (— I 2^^3 — <T$rg), 

A2 = ^9^5^18 • • • ( — '^1W< 
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Substituting m'Q^i for g^f and dividing by ni'Dja, we have 

»t'-'cg^ = - (i - A^)Qz + ^Ge - f ^a, 

m>-'ag,= eQ,-{6^ + i<T)Q„ (17) 

and, for succeeding terms, 

^3i = 4-'<23. = Y^T^ '^'^^si- (18) 

6. For the sine terms not involving multiples of 3^, we may assume /i = ^ 
and ;3 = I ; whence, substituting in (10) and separating by intederminate coeffi- 
cients, 

^'=7~7i*»'Qi- (19) 

y ' 

7. Instead of (3), we may write 

ndt ' ^ >^ ^a n' — n 

= 1/(1 — v) + ni'I,{Si cos «<^ + Ti sin ih), (20) 

in which St and Zj are assumed to be functions of :y only. We have, also, 

If we subtract (7) from (20), neglecting powers and products of bs, dv = v — e^, 
and 

» dw dw dw, 

ndt ndt ^dd 
we obtain 

2s, ^Ss + s,'d ^^= - 1 (I + ^ ,^) ^v + m'l\{S, cos id + T, sin id), (21) 

in which 

s, -~ =1 —\e^ — ^ cos 3^ -I- ^ cos 66, 

si^ = 1 -\- ^e^ -\- 26 cos 3^ — J^ cos 6^. 
Assuming 



^ ~-j = ^i{»i cos ?'(? + hi sin «^), (22) 
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in which «< and hf are functions oi'/j only, substituting in (21), and separating by- 
indeterminate coefficients, we have, approximately, 

^«, + (i + f '^)«o H- ^«3 — l<^n^ = — Hi + J^)'^". 
— f«e + 2£'«„ + «3 + ^«8 = w'^a, 

the solution of which gives 

«„ = — J(i + ^)^i' — f ^a, — ^w'53 + \^m'S^, 

«3= ^^1'+ 3^«6+ m'S^— em' St, (23) 

«« = — 2^8 + fn'S^ ; 

or, substituting from equations (12), 

,«'-•«„=- Hi - J^)i'„ + IrP, - eS, + f^S^, 

m>-'-n^= eP^— eP^^ S, — eS,, (23') 

Also, 

2(1 - 0^3 - ^.fe + (I + ■|-^)/'3 + ^^^6 = ^t'T„ 

— egi + 2^8 + eh^ -\- h^ =»i'T^, (24) 

2^9 — i^'^3 + ^^6 + -^9 = =^'^9 ; 

whence 

hi = w'T; — 2(1 — 1^^3 + 3<f8, 

K = m'T^+ 3<r3— 2^8, (24') 

^g = tn' 7; — I- r'^j + 2^^, — 2^9 ; 

or, substituting from equations (17), 

m'-'ah, = ^7; + 2(1 - J/.^) Q, - 2eQ, + 3^^,, 

For the remaining terms, we may write, in general, 

2a^ — ea,+i — eai_3 + «i + ^«i+3 + ^«i-3 = WSi, (25) 

2gi — ^£'i+3 — ^^i-3 + hi + eki+3 + eAi_, = m'Ti, (26) 



52 STONE. ON THE MOTION OF HYPERION. 

in which 

<§"•-.!= g-i + i, «._3 = — »_i + s; 

whence, approximately, 

Wi±3 = ^'Si±i — 2ai+^, 
hi±-i — m'Ti+i — 2gi+i ; 
which, substituted in (25) and (26), give 

«i = m'Si — cm'Si+:i — em'Si_i — 2a, + s^i+j + 3f«i_3, (27) 

hi = wTi — eni'Ti^^ — em'Ti_i — 2^; + 3^^i+3 + 3^.ri-3, (28) 

in which 

8. If we omit periodic terms, we may write 



{dw'\ __ (dwj^ { A^_\ 

\^ndt} l3«'<^J I ndt] 



or ? = I 4- ^»«', (29) 

in which A is the constant part oi m'~^n^, which is given by the first of equa- 
tions (23'). We may also put 

a = gdii = 9 Kz_^' _ , = Bm', (30) 

in which B is obtained from the second of equations (12). 
Putting 

l{n'-n ) ^_zn'-~4n 

n n ^'' ' 

and substituting in (30), we have 

(i + ff{\ + AmJ= I + Bm'; 

or, with sufficient accuracy, 

2?" + r^ + 2Am.' = Bm'; 
whence 



STONE. ON THE MOTION OF HYPERION. S3 

9. We may now put 

ndt n *" y^ndtj » 

= 5 + "'+[f^].+(*S]/ (33) 



in which 



n,, ^ n,, — . Am'. 



The subscript indicates that only terms involving that quantity are considered. 
Multiplying by ndi and integrating, 

w = Wo + «/ + -, fn^ri + -^ («/,)» + -J^— fid^-'^] dd; (34) 

4« — 3«'J ^ ' T,{n' — n)^'' n' ~ nj y, ndt) e ^^^' 

or, with sufficient accuracy, 

-^nt- ^Jn,dri + (i-Am'-r) («/,), + s/ [ ^ ^] ^ ^^. (34') 



Since 



where 



g — 4^ — 3«' ^ _ 3r 



we have, by integration, 

'? = g(?H-g„; 
whence 

cosy:y cos id = J cos {id +jyj) + J cos {id —jrj) 

= i cos i{i + jg) e +7gj + J cos [(/ -yg) d -yg„] , 

and 

/ cos jrj cos «^ af(? = J . ,' . sin (/(? + jrj) + A ,— L^. sin {id —jrj) 
J t -\-jg " t ~jg ^ -' '^ 

= ^r^ ^^o^^"? «'" ^'^ - ,-. j^-2g2 siny^? cos ^^. (35) 
Similarly 

J sin jyj cos ?(? (/^ = .^ _^ sin jrj sin id + ^.^ j^^.^ ^ cos 7;; cos id, (36) 
J cosy;y sin id dd = — j^ ^ cos jrj cos id — ^., ^.^ ^ sin jrj sin «"^, (37) 
J sin /if sin id dd = — ^., _' ., . sin 7;; cos id + .^ jj^.^^^ cos jrj sin «^. (38) 
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Since, however, g is small and is multiplied by small coefficients, it may, in 
general, be neglected ; whence (34') becomes 



w 



I C * I 

= Wd + «/ — - I n^rj-{- (i — Am' — r)(^/)» + 3-^V- («,sin zW — i^,- cos /^) . (34") 



Similarly, in integrating Srdd, the coefficients of sines and cosines of multiples 
of 6 may, in general, be considered as constant. 
10. Differentiating 

V{i - V) = ^(i - V.) + ^ ^4^/(^^)o'^^ (39) 

with regard to ly, we have 

1 /t 1 i,\^^ ^z. I i^N^^ I d^^' in' n' .-^ 

d-q dri ' dfj ya 4« — 3«' ^ •^° 



or 



W. = -^fli^'(M.+K)| + if], (40) 

and may assume 

e=e^{i +^iCos)y +/isin^ + ^2Cos2^ +/jsin2^ + ...). (40 

in which e„, ^j./i, etc. are constants. 

11. The coefficients of sines and cosines of multiples of yj in e, 8v, at, tii, gi, 
and hi are found by substituting (41) in equations (12), (17), and (23)-(26), and 
separating by indeterminate coefficients. With assumed values of e^, /J,, P^, etc., 
the numerical values of B and of the coefficients in the expression for ej e^ are 
first obtained from the second of equations (12); the numerical values of the 
other coefficients are then found by substituting these in the remainder of the 
equations mentioned. 

12. For the computation of the forces we have the well-known formulae, 

p2 __ ^/2 gj^2 (^^1 — w) -\- [r' cos [w' — w) — r]^, 
t I I 



p^ r'^ ' 



R = hr' cos {-w' — w) — -^ , 
5' = hr' sin («/' — w) , 

in which p is the distance between the two satellites, and r' and w' are the radius 
vector and longitude in orbit of Titan. 

Assuming an approximate orbit, a^R and Sr are computed for different values 
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of fj and d. In so doing ij and 6 may be treated as constants. The assumption 
however, that ly is constant is equivalent to saying that the conjunction point is 
stationary, or, in other words, that the mean motions of the two satellites are 
commensurable. If 4 be the mean time of conjunction of the two satellites, we 
have 

^ = /'— /— /i=I' — 1— yi, 
in which 

\-^ri^l{n' —n)(t~t,) 
and 

l'=^+4(«' — «)(/ — 4) 

are the mean longitudes of Hyperion and Titan on the assumption that ^ is con- 
stant. 

Having determined the values of c?R and Sr for individual values of G and 
ly, the former are expressed in series of sines and cosines of the latter by means 
of mechanical quadratures. With the values of c^R and Sr thus expressed /"(, 
Qi, Si, 7^, are obtained and substituted in equations (12), (17), and (23) -(26), 
whence are derived the inequalities in r and w. The problem is solved when we 
have found expressions for r and w which agree with the assumed values of those 
quantities. The value of {Sr\ found by mechanical quadratures must also agree 
with that obtained by means of equation (40). 

13. Assuming the perisaturnium of Titan as the origin of longitudes; putting 

X = o, 

«= 16°. 9199, «'^ 22°. 5770, 

^ = ^" = 0.01, e' ^ 0.0284, 

log a ^ 0.08351, a'^i; 

and neglecting ds and d —j , the following values of c^R and Sr were obtained 
for ly ^ o and ly = 1 80° : — 



S6 
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'? = 


= 0° 




"/ = 


180° 


f) 






























a'R 


Sr 




a'R 


Sr 




CO 


— 


12.64 


0.00 


— 


17.20 


0.00 


7-5 


— 


8.17 


+ 3-23 


— 


1 1.22 


+ 4-7' 


15.0 


— 


4-32 


+ 2.02 


— 


576 


+ 342 


22.5 


— 


2.83. 


+ 1.07 


— 


361 


+ 1 87 


30.0 


— 


215 


+ 0.52 


— 


2.68 


+ 091 


37-5 


— 


1.81 


+ o,J4 


— 


2.17 


+ 035 


45.0 


— 


■63 


+ 0.13 


— 


1.S4 


-\~ 0.02 


52.5 


— 


1.56 


+ 0.07 


— 


..63 


— 0.12 


60.0 


— 


1.58 


+ 0.0 1 


— 


1.42 


— 0.25 


67.5 


— 


1.52 


— 0.1 1 


— 


1.28 


-0.33 


75.0 


— 


1.40 


— 0.33 


— 


1.14 


— 0.44 


82.5 


— 


1. 18 


-058 


— 


0.98 


— 0.54 


90.0 


— 


0.86 


— 0.85 


— 


0.78 


— 0.72 


97-5 


— 


0.46 


— 1.09 


— 


0.57 


--0.81 


105.0 


— 


0.12 


- 1.09 


— 


0.28 


— 1. 00 


112.5 


+ 


0.20 


— 1.06 


+ 


003 


- i.oS 


120.0 


+ 


0.46 


— 96 


+ 


0.34 


— 1.06 


127.5 


+ 


0.66 


— 0.84 


+ 


0.67 


— 0.96 


135-0 


+ 


0.81 


— 0.68 


+ 


0.93 


— 0.77 


142.5 


+ 


0.93 


— 0.53 


+ 


1.09 


— 0.53 


150.0 


+ 


1.03 


— 0.39 


+ 


1. 15 


— 0.31 


>57-5 


+ 


I. II 


— 0.27 


+ 


1. 14 


— 0.15 


165.0 


+ 


1. 16 


— 0.17 


+ 


I. II 


— 0.06 


172.5 


+ 


1. 19 


— 0.08 


+ 


1.07 


— 0.02 


180.0 


+ 


1.06 


0.00 


-1- 


1. 21 


0.00 



For negative values of d, the signs of Sr change, but not those of a^R. 
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Considering terms involving multiples of 3</ only, the following values of 
a^R and Sr give 

a^R^ — 1.326— 1.224 cos S'^ — 0.847 cos 6(^ — 0.597 cos 9^ — 0.323 cos 12^ 

+(0.167 + 0.276 cos 3^^ + 0.158 cos 6^ + 0.084 cos 9^ + 0.043 cos i2^)cos;y 

Sr= + 0.633 sin 3^ +0.534 sin 6Q +0.356 sin 9^ + 0.230 sin 12^ 

+(0.103 sin 3^ — 0.091 sin 66 — 0.069 sin 9^ — 0.039 ^'^i i2^)cos)f. 

We have, also, 

1/(1 — v) n' 

-!:--i ^ —, = 1.202 

3j/« «' — n 

and 

- = -J ^ 1.015 — 0.303 cos 3^ -I- 0.045 cos 6^ — 0.007 cos 9^ + o.ooi cos \2d\ 

o S f 

whence 

P„ = — 1. 108 + 0.132 cos -fi, 

/"s — — 2 709 + 0.5 16 cos //, 

F^ = — 1.225 + 0.221 cos -tj, 

Pj = — 0.806 + 0.1 30 cos if, 

Substituting in the first of equations (23'), we have 

vi'' ' «„ =: 0.639 — 0.069 cos ij ; 
whence 

A = 0.639 and m'— ' «, = — 0.069 cos fj. 

The assumed values of n and «' give i2fJ/.< = 0.0081 ; whence, by the second 

of equations (12), 

gm'-^dn . e = 2.883 — 0.526 cos q ; 

whence, assuming ^ = ^ (i -^ e^ cos •//), we have 

^ = 9w'-'o>= 28.83 



s,= 


— 0.765 + 0.124 COS )f. 


s,= 


— 0.322 + 0.054 cos ■>}, 


s,^ 


— 0.145 + 0.029 cos :y. 



and 
i. e. 

We have, also, 



e^= —0.183; 
f = o. 1 000 — 0.0 1 83 cos fl. 

in' — 4« = o°.05i4, 
;- = 0.003038 ; 
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whence, substituting in (32), 

, I 

4528' 

and the long period inequality in (34") becomes 

I C J 0.o6q X 57°.3 sin rt o o • 
— — I nav = — ;:^2 iZ — 2 ^ = o°.287 sin yj. 

rJ ' ' 0.003038 X 4528 ' ' 

The terms in d, 28, 46, and $6 give 
a^R— — 3.23 cos 6 — 1. 31 cos 2d — 1.31 cos4<^ — 1.98 cos 5</ 

-^ (0.38 cos ^ + 0.31 COS 2O + 0.21 cos 4.0 + 0.07 COS 5</) cos Tj, 
Sr = — 0.44 sin + 0.76 sin 26 + 0.39 sin 4^ + 0.57 sin 5^ 

+( — 0.06 sin — 0.06 sin 2O — 0.14 sin 46 — 0.17 sin 5^) cos vj; 
whence 

-^1 = + 0.46 + 0.39 cos if, 5i = + 1.59 -f 0.22 cos Yj, 

P^= — 4.46 + o. 1 3 cos -/f , ^2 = — 1 .38 + o. 1 1 cos --f , 

/*4 = — 2.47 + o. I s cos Tj, 5^ = — 0.35 -\- o. 1 2 cos if ^ 

/'s = — 2.29 -}- o.i9cos;y; 5^ = —0.41 +0.12 cos ;f . 

The resulting values of a^ and «, are 
^1 — — 0.0002 + 0.000 1 cos -q, «, = -|- 0.0004 — 0.0002 cos 7), 

a.2= — 0.00 1 8 + 0.000 1 cos -q, «2 =" + 0.0034 — 0.000 1 cos -y, 

«;, = 0.0000, «3= — 0.0002, 

«4 =^ + 0.0005 — 0.000 1 cos yj, Mi = — 0.0012 + 0.0002 cos /j, 

«6 = + 0.0006, «5 = — 0.0018 + o.oooi cos :y, 

«6 = ^^ 0.000 1 , «j = — 0.0003 ; 

^7=: — 0.0001 ; 
whence, neglecting terms involving^; and those containing both rj and d, 
r =: a(s, + OS) = as, — o.OOOi cos Q — 0.0022 cos 26 

+ 0.0006 cos 4^ -f 0.0007 cos S^' + 0.0001 cos 6# 
— 0.0001 cos yd. 
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We have, also, 

«/, = w„ + nt — {2e — ■J-^)sin 3<^ -\- \^ s\n 66 — -J-I^^sin 9^ + -^^^sin 126, 

in which 3^ is the mean anomaly counted from a/i^saturnium, and 

{Am' + y) {w,)e = — 2e (Am' + y) sin 3^ 

=: ( — 0.0006 + o.oooi cos -f) sin 3^; 

whence, neglecting terms involving ^j and those containing both rj and 0, equation 
(34") gives 

zi> =z w, -\- o°.29 sin r^ -\- o°.i7 sin f) + o°.29 sin 2t) — o°.oi sin 3^ 

— o°.05 sin 4.0 — o°.o6 sin $d — o°.oi sin 60. 

14. For )j ^ 90°, the forces were computed on the same hypothesis as 
before, and also on the assumption that w contains the additional term, 

Jw = — A3 cos 3^ = — 2°. 50 cos 3^. 

The corresponding term in r is, approximately, 

Jr = ag^ sin 3^ ^ — \cilh ^'n 3^ = — 0.0264 sin -xft. 

For fj = 270°, when / ^ o, the numerical values of the forces are the same 
as for /j = 90°; except that the sine terms in a^R and the cosine terms in Sr 
change sign. 

The values of a^R and Sr obtained for r^ = 90° are as follows: — 



6o 
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Jw ■ 


= 


Jw:= 2 


°.5 cos id 


e 




















a'R 


Sr 


a-'R 


Sr 





— 14.08 


+ 1-87 


— 12.81 


+ 2.93 


2-5 


- 12.31 


+ 3-47 


— 10.77 


+ 3-94 


5 


— 10.03 


+ 4-1 1 


— 8.69 


-h 4-14 


7-5 


- 8.03 


+ 4-05 


— 6.98 


+ 3-86 


lO 


— 644 


+ 3-65 


- 5-69 


+ 341 


15 


— 4-39 


+ 2.67 


— 4-07 


+ 2.50 


20 


— 3-34 


+ 1.86 


— 319 


+ 1.80 


25 


- 2.75 


+ 1.29 


— 2.69 


+ 1.31 


30 


- 2.38 


+ 0.89 


- 2.38 


+ 0.97 


35 


- 2.15 


+ 0.62 


- 2.17 


+ 0.72 


40 


— 1.98 


+ 0.41 


— 2.04 


+ 0.54 


45 


— 1.88 


+ 0.25 


- 195 


+ 0,39 


50 


- 179 


+ 0.13 


- 1.87 


-h 0.26 


55 


- 1.72 


+ 0.02 


— 1.81 


+ 0.14 


60 


- 1.63 


— O.IO 


- 1-73 


— 0.0 1 


70 


— 1.41 


-0.35 


— 1.46 


— 0.30 


80 


— 1.09 


— 061 


— 113 


— 0.62 


90 


— 0.71 


-0.83 


— 0.70 


— 0.86 


100 


': — 0.33 


— 0.96 


— 0.29 


— I.OI 


no 


+ 0.02 


— I.OI 


+ 0.08 


— 1.04 


120 


+ 0.35 


— 0.98 


+ 0.41 


— 0.95 


130 


+ 0.67 


— 0.87 


+ 0.72 


— 0.82 


140 


+ 0.94 


— 0.65 


f 0.96 


— 0.62 


150 


+ I-I4 


— 0.40 


+ 1.13 


— 0.40 


160 


+ 1.20 


— 0.19 


+ 1.20 


— 0.19 


170 


+ 1. 19 


— 0.02 


+ 1. 18 


-- 0.06 
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e 


Aw -- 


= 


Aw = — 2 


°.S cos 3/? 


a^R 


Sr 


a^R 


Sr 


1 80 


+ 1.13 


+ o.os 


+ I-13 


+ 0.07 


190 


+ 1.09 


+ O.IO 


+ 1.09 


+ 0.0 1 


200 


+ 1-07 


+ 0.18 


+ 1-07 


+ o.is 


210 


+ 1.04 


+ 0.31 


+ 1-05 


+ 0.32 


220 


+ 0.98 


+ 0.55 


+ 0.97 


+ 0.56 


230 


+ 0.82 


+ 0.81 


+ 0.78 


+ 0.85 


240 


+ 0.47 


+ 1.06 


+ 0.42 


+ 1.08 


250 


0.00 


+ i-iS 


— 0.06 


+ 1-13 


260 


— 0.50 


+ 1. 01 


- o.ss 


+ 0.97 


270 


— 0.90 


+ 0.73 


— 0.93 


+ 0.69 


280 


- 1. 19 


+ 0-41 


- 1-15 


+ 0.41 


290 


- 1.30 


+ 0.21 


— 1-23 


-+ 0.26 


300 


- 1.36 


+ 0.13 


— 1.28 


+ 0.23 


305 


— 1.40 


+ 0.12 


— 1-32 


+ 0.22 


310 


— 1.48 


-f O.II 


— 1.40 


+ 0.21 


31S 


— 1.60 


+ 0.07 


- 1-53 


+ 0.17 


320 


— 1-77 


— 0.02 


— 1-71 


+ 0.08 


32s 


— 2.04 


— 0.20 


— 1-99 


— O.II 


330 


— 2.40 


— 0.50 


- 2.39 


— 0.43 


335 


— 2.94 


— 0.98 


— 2.98 


— 0.94 


340 


- 3-83 


— 1.70 


- 3-98 


— 1.67 


345 


— 5-49 


— 2.70 


— 6.02 


-2.71 


350 


— 8.65 


-3-52 


— 9-44 


— 3-25 


352.5 


— 10.89 


— 3-35 


— 11.71 


— 2.66 


355 


— 13-12 


-2.31 


- 13-54 


— 1. 12 


357-5 


— 14-43 


-0.35 


- 13-98 


+ 1.04 
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This table was completed by interpolation so as to contain values of the 
forces for every two and a half degrees of d ; whence, considering terms involving 
multiples of 3^ only, we have for dw = 0, dr = 0, 

a^R= — 1.287 — i.i68cos3^ — 0.8oicos6^ — 0.489 cos 9^ — 0.294 cos 12^ 

+(0.059 sin 3^ +0.134 sin 6^ + 0.131 sin 9^ +0.109 ^in I2^)sinj;, 

Sr^= 0.654 sin id + 0.506 sin 60 + 0.408 sin gd + 0.214 sin 12^ 

+ 0.127 sin 15^ + 0.072 sin 18^ 

4 (0.077 + 0.091 cos 3^ + o.o89cos6^ + 0.094 cos 9^ + 0.078 cos 12^ 

+ 0.064 cos 15^ +0.047 cos 18^) sin;y; 
whence 

Q^ = 0.261 sin Yj, 

Q^ = 0.200 sin -f] , 

Q^ = 0.186 sin ;;, 

012=0.141 sin -fj. 

T^, Zj, etc., are negligible. 

In the same manner, for dw = — 2°. 5 cos 7,6 , dr = — 0.0264 sin 3^ , 

d^R^^ — I 259 — I. II 3 cos 3^ — 0.743 cos 6^ — 0.429 cos 9^ — 0.250 cos 12^ 

+(0. no sin 2,d +0.227 sin 6^ +0.224 sin 9^ +0.1 80 sin 12^) sin jy, 

Sr== 0.635 sin 3^ + 0.468 sin 6^ +0.372 sin 9^ + 0.176 sin 126 

+ 0.093 sin 1 5^ + 0.045 sin 18^ 

+(0.128 +0.126 cos id +o.i58cos6^ + o.i59cos9^ +0.133 cos 12^ 

+ 0.100 cos 15^ + 0.074 cos 1 8^) sin :^, 

^^ ' —, U= — 0.763 cos 3^ — 0.28 1 cos 6^ — 0.1 26 cos 9^ — 0.05 3 cos 12^ 

x/a n' — n /jo y jo 

— 0.023 cos 15^ — 0.009 cos 18^ 

+(0.151 sin 3^ +0.095 sin 6^ + 0.064 .sin 9^ + 0.040 sin 12^ 
+ 0.024 sin 15^ +0.015 sin 18^) sin ly. 
Multiplying the latter by 
2s~^ = 2.030 — 0.606 cos 3^ + 0.090 cos 6d — 0.014 cos 9^ + 0.002 cos 12^ 
+(0.158 sin 2)6 — 0.008 sin 60) sin tj, 
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we have 

U' = (0.301 sin 3^ + 0.077 ■'*'" 6^ + 0.085 sin 9^ + 0.052 sin \2ff) sin )j; 

and, adding c?R, 

03 = 0.411 sin ;y, 

^6 = 0.304 sin jy, 

(2g = 0.309 sin )y, 

(2i2= 0.232 sin -q. 

In the present case, in a first approximation, the values of Q.^, Q^, etc., may 
be considered as depending on the value of Jzc for = o. 

The above values give, by means of (24') and (34"), for r^ = 90°, 

Assumed. Computed. 

Jw= — 2°.s cos 3^, Jw= +o°.287 — i°.47i cos 3^ + 0°. 195 cos 6^ 

— 0° .024 cos gd + o°.oo2 cos 12^ , 
^ — I °.Oi I , for ^ := o, 

Jw = 0, Jw= +o°.287 — o°.929cos 3^ + o°.i23 cos 6^ 

— o°.oi4cos gd +o°.oo2 cos 12^, 

= — 0°.53i, for d = o; 

whence, by interpolation, Jw (assumed) = Jw(computed) = — o°.657, for ^ = o. 
Hence, approximately, 

03 =z 0.300 sin 7], 

08 ^0.227 sin 7], 

09 = 0.218 sin V], 
0,2=0.165 sin r;, 

Jw — 0° .2g — 1° .oy cos 3^ + 0°. 14 cos6(?— o°.02 cosgd, 
Jr^=: — o.oioo sin 3^ + o.ooio sin 66 — o.oooi sin gd . 

For 7} = 270° and .^ = O, the values of these inequalities are the same as 
those for !y =^ 90°, but with opposite signs. 

In a similar manner, a^R, Sr, Aw, and Ar may be found for other values 
of ly, and the coefficients of the sines and cosines of multiples of d expressed as 
periodic functions of that quantity (vj) by means of mechanical quadratures. 
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15. With the values of the co-ordinates thus found more accurate values 
may be found by a repetition of the process described in the preceding pages ; 
but the solution cannot be considered as even approaching completion until the 
value of X is obtained by the comparison of an assumed orbit with observations 
made near conjunction, and distributed, if possible, over the whole of the nineteen 
years during which vj makes a cycle of 360°. 

It may be noted that (Sr\ as given by equation (40) is nearly proportional 
to de^jdrj, and, therefore, is liable to converge slowly, since e^ converges more 
slowly than e, and de^jdrj converges more slowly than e^. Hence, the values of 
{Sr\ obtained for 90° and 270° are not sufficient to give, approximately even, 

the coefficient of sin rj in that quantity. On the other hand, — — I «,</jj con- 
verges more rapidly than «, ; whence it is probable that the value of that in- 
equality may be derived with considerable accuracy from the forces obtained for 
;y = O and )y =: 180°; but, in computing those forces the inequality in e should 
not be neglected, except in a first approximation. 

16. The method here suggested may be readily extended to cases in which 
the disturbed and disturbing bodies do not move in coplanar orbits. The a.ste- 
roids furnish quite a number of interesting cases of mean motions nearly com- 
mensurate with that of Jupiter. As examples maybe mentioned (153) Hilda 
and (190) Ismene, whose mean motions are each approximately three halves that 
of Jupiter. 



